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Abstract 


27 The chemical link between isoprene and formaldehyde (HCHO) is a strong, non-linear function of NO x (= 

28 NO + N0 2 ). This relationship is a linchpin for top-down isoprene emission inventory verification from 

29 orbital HCHO column observations. It is also a benchmark for overall mechanism performance with 

30 regard to VOC oxidation. Using a comprehensive suite of airborne in situ observations over the 

31 Southeast U.S., we quantify HCHO production across the urban-rural spectrum. Analysis of isoprene and 

32 its major first-generation oxidation products allows us to define both a "prompt" yield of HCHO 

33 (molecules of HCHO produced per molecule of freshly-emitted isoprene) and the background HCHO 

34 mixing ratio (from oxidation of longer-lived hydrocarbons). Over the range of observed NO x values 

35 (roughly 0.1 - 2 ppbv), the prompt yield increases by a factor of 3 (from 0.3 to 0.9), while background 

36 HCHO increases by more than a factor of 2 (from 1.5 to 3.3 ppbv). We apply the same method to 

37 evaluate the performance of both a global chemical transport model (AM3) and a measurement- 

38 constrained 0-D chemical box model. Both models reproduce the NO x dependence of the prompt HCHO 

39 yield, illustrating that models with updated isoprene oxidation mechanisms can adequately capture the 

40 link between HCHO and recent isoprene emissions. On the other hand, both models under-estimate 

41 background HCHO mixing ratios, suggesting missing HCHO precursors, inadequate representation of 

42 later-generation isoprene degradation and/or under-estimated hydroxyl radical concentrations. 

43 Moreover, we find that the total organic peroxy radical production rate is essentially independent of 

44 NO x , as the increase in oxidizing capacity with NO x is largely balanced by a decrease in VOC reactivity. 

45 Thus, the observed NO x dependence of HCHO mainly reflects the changing fate of organic peroxy 

46 radicals. 

47 

48 1 . Introduction 

49 Formaldehyde (HCHO) is a ubiquitous byproduct of volatile organic compound (VOC) oxidation. While 

50 methane is the principal HCHO precursor in remote regions, larger VOC are the main source over 

51 continents. HCHO is also directly emitted via biomass burning (Lee et al., 1997), fossil fuel combustion 

52 (Luecken et al., 2012), natural gas flaring (Knighton et al., 2012), ethanol refining (de Gouw et al., 2015), 

53 possibly vegetation (DiGangi et al., 2011) and agricultural activity (Kaiser et al., 2015a), but chemical 

54 production dominates the global budget (Fortems-Cheiney et al., 2012). Photolysis and reaction with OH 
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55 destroy HCHO with a characteristic lifetime of several hours during midday, implying that the HCHO 

56 abundance reflects recent hydrocarbon oxidation. 

57 Globally, isoprene is the main precursor of near-surface HCHO. A highly reactive diene emitted 

58 by vegetation, isoprene comprises roughly one third of all non-methane VOC emissions (Guenther et al., 

59 2012). Oxidation of isoprene in the presence of nitrogen oxides (NO x = NO + N0 2 ) stimulates the 

60 production of ozone (Trainer et al., 1987) and organic aerosol precursors (Xu et al., 2015), impacting air 

61 quality and climate in many continental regions. Despite the central role of isoprene, biogenic emission 

62 inventories struggle to accurately represent the spatiotemporal variability of isoprene emissions, with 

63 model-measurement discrepancies and differences among emission inventories approaching a factor of 

64 2 or more (Carlton and Baker, 2011; Warneke et al., 2010). Such differences directly impact predicted 

65 ozone and aerosol distributions (Hogrefe et al., 2011). 

66 Numerous studies have applied satellite-based HCHO column observations as a top-down 

67 constraint on isoprene emissions (see Kefauver et al. (2014) for a review). Typically, a chemical transport 

68 model is employed to relate HCHO column concentrations to isoprene emission strength. Early studies 

69 utilized linear steady-state relationships (Palmer et al., 2003), while recent computational advances have 

70 permitted full inversions that more fully account for transport, multiple sources and varying chemical 

71 regimes (Fortems-Cheiney et al., 2012). Such techniques have informed isoprene emission inventories in 

72 North America (Abbot et al., 2003; Millet et al., 2008; Millet et al., 2006; Palmer et al., 2006; Palmer et 

73 al., 2003), South America (Barkley et al., 2013; Barkley et al., 2008), Europe (Curci et al., 2010; Dufour et 

74 al., 2009), Africa (Marais et al., 2012), Asia (Fu et al., 2007; Stavrakou et al., 2014), and globally 

75 (Fortems-Cheiney et al., 2012; Shim et al., 2005; Stavrakou et al., 2009). Future geostationary 

76 observations, such as the NASA Tropospheric Emissions: Monitoring of Pollution (TEMPO, 

77 http://science.nasa.gov/missions/tempo/) mission, will permit an even more detailed investigation of 

78 the spatial and temporal variability of isoprene emissions and other VOC sources. 

79 Chemistry dictates the relationship between HCHO columns and underlying isoprene emissions. 

80 Many of the above-listed studies apply 0-D box model calculations to evaluate the yield of HCHO from 

81 isoprene as a function of oxidation time, NO x regime and chemical mechanism. In all cases, it is found 

82 that NO x enhances both the production rate and ultimate yield of HCHO. Slower production at lower 

83 NO x can lead to "smearing," whereby HCHO production is displaced relative to the isoprene source. 

84 Palmer et al. (2003) define a characteristic smearing length scale, which can range from 10 to 100 km or 

85 more. Furthermore, accumulation of oxygenated VOC over multiple generations of isoprene 
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86 degradation can contribute to substantial background HCHO production, which is not directly linked 

87 with fresh isoprene emissions. Long-lived primary anthropogenic or biogenic emissions, like methane 

88 and methanol, can also contribute to this background. Background column concentrations are typically 

89 on the order of 5 x 10 15 cm 2 , equating to 20% or more of the isoprene-driven HCHO column 

90 enhancement (Barkley et al., 2013; Millet et al., 2006). A wave of recent theoretical (Peeters et al., 2014; 

91 Peeters and Muller, 2010; Peeters et al., 1999), laboratory (Crounse et al., 2012; Crounse et al., 2011; 

92 Paulot et al., 2009a; Paulot et al., 2009b) and field (Mao et al., 2012) research has highlighted 

93 shortcomings in low-NO x isoprene oxidation schemes. Such issues translate directly into top-down 

94 emission estimates; for example, Marais et al. (2012) report an uncertainty of 40% in OMI-derived 

95 African isoprene emissions at high-NO x and 40-90% at low-NO x . Coarse resolution of averaged satellite 

96 observations and model simulations (typically 1° x 1° or more) has partly mitigated these problems in 

97 prior work, as variability in NO x -dependent smearing and background production is averaged out. A 

98 more careful treatment will be needed to harness the enhanced resolution of near-future orbital 

99 observations (e.g., 8 x 4.5 km 2 for TEMPO), especially since these measurements will include diurnal 

100 variability. 

101 Here, we use a comprehensive set of in situ observations to quantify the impact of NO x on the 

102 isoprene-HCHO chemical link. Using isoprene and its unique first-generation products, we segregate 

103 HCHO into two categories. The first, defined as "prompt" HCHO, is produced from fresh isoprene 

104 emissions (on a timescale of less than a day) and retains the signature of isoprene emission source 

105 strength. The second category is "background" HCHO stemming from oxidation of longer-lived isoprene 

106 oxidation products and other VOC. We examine the NO x dependence of both quantities. Applying the 

107 same method to 0-D and global model simulations, we evaluate the ability of current chemical 

108 mechanisms to replicate the observed trends. Box model results are also used to elucidate the 

109 mechanistic underpinnings of the NO x influence on HCHO production. 

110 

111 2. SENEX Observations 

112 The Southeast Nexus (SENEX) mission was an airborne campaign designed to examine the interaction of 

113 natural and anthropogenic emissions. During June and July of 2013, the NOAA WP-3D aircraft logged 

114 114 flight hours over 18 research flights in a range of environments throughout the southeast United 

115 States, including urban centers, power plant plumes, natural gas extraction regions, agricultural areas 
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116 and forests. The payload included a suite of gas- and particle-phase instrumentation (Warneke et al., in 

117 preparation, 2015); details and data are accessible on the SENEX website 

118 (http://www.esrl.noaa.gov/csd/projects/senex/). Here we utilize observations of HCHO, isoprene, 

119 methyl vinyl ketone (MVK), methacrolein (MACR), NO and N0 2 . HCHO was measured at 1 Hz by the 

120 NASA In Situ Airborne Formaldehyde (ISAF) instrument, which relies on the laser-induced fluorescence 

121 technique and has an accuracy of ±10% (Cazorla et al., 2015). Isoprene, MVK and MACR were measured 

122 by both a quadrupole proton transfer reaction mass spectrometer (PTR-MS) and the NOAA improved 

123 whole-air sampler (iWAS) with offline gas chromatography. The PTR-MS (de Gouw and Warneke, 2007) 

124 has a stated accuracy of 20% and sequentially sampled masses for isoprene (m/z +69) and the sum of 

125 MVK and MACR (m/z +71) for 1 s each with a duty cycle of 14 s. The iWAS measurement uncertainty for 

126 speciated MVK and MACR is 20% (de Gouw et al., 2015). NO and N0 2 were measured at 1 Hz via 

127 chemiluminescence coupled with a photolytic N0 2 converter (Pollack et al., 2010; Ryerson et al., 1999) 

128 with an accuracy of 5%. Data are filtered to include only daytime boundary layer conditions (solar zenith 

129 angle < 60°, radar altitude < 1 km). Influence from biomass burning (acetonitrile > 210 pptv and CO > 

130 300 ppbv) and very fresh power plant plumes (log(NO x ) values exceeding a mean + 3a threshold) are 

131 also removed. This procedure excludes 50% of all fast (1 Hz) data. After accounting for missing data, we 

132 retain 8435 1 Hz data points and 81 iWAS samples. 

133 Measurements of MVK and MACR can include a positive bias from conversion of isoprene 

134 hydroxyhydroperoxides (ISOPOOH) on hot metal surfaces in the sampling system (Liu et al., 2013; 

135 Rivera-Rios et al., 2014). Theoretically, this mechanism could give rise to an analogous artifact in HCHO 

136 observations. ISOPOOH mixing ratios of roughly 0 to 2 ppbv were observed during SENEX (see 

137 supporting information (SI)). It is difficult to quantify the magnitude of any such interference from field 

138 observations alone. Based on a comparison to observations of other isoprene oxidation products and to 

139 0-D box model results (SI), we argue that such artifacts are negligibly small in the PTR-MS and ISAF 

140 observations for SENEX. We cannot rule out a potential positive bias in the iWAS MVK measurement; 

141 nonetheless, as we show below, the correspondence between observed MVK and MACR mixing ratios is 

142 consistent with our current understanding of isoprene oxidation. 

143 SENEX sampled a wide spectrum of chemical regimes (Figure 1). For the daytime boundary-layer 

144 observations presented here, maximum 1 Hz isoprene and NO mixing ratios respectively reach 8.1 and 

145 95 ppbv, while minima are less than a few pptv. The distributions of both isoprene and NO observations 

146 are approximately log-normal (top and right panels of Fig. 1), peaking at 1.5 ppbv and 50 pptv, 
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147 respectively. Though these distributions may be biased towards areas of urban influence, the range of 

148 environments encountered during SENEX is representative of the Southeast U.S. summertime boundary 

149 layer. The long tail at the low end of the isoprene distribution is mostly associated with regions lacking 

150 significant tree cover, notably Illinois and Indiana, where isoprene emissions are generally lower. The NO 

151 distribution spans two orders of magnitude (10 - 1000 pptv), over which radical chemistry changes 

152 markedly. At NO mixing ratios of a few hundred pptv or more, organic peroxy radicals (R0 2 ) react mostly 

153 with NO. At low NO (10's of pptv or less), reaction with H0 2 , other R0 2 and isomerization dominate. The 

154 bulk of the NO x distribution lies in a transition region for radical chemistry, making this dataset ideal for 

155 probing the anthropogenic influence on biogenic VOC oxidation. 

156 HCHO mixing ratios (color shading in Fig. 1) range from 0.8 to 14 ppbv with a mean value of 4.3 

157 ppbv. HCHO is most abundant in regions where both isoprene and NO x are elevated. High NO x is often 

158 accompanied by increased concentrations of anthropogenic VOC; however, constrained box-model 

159 calculations demonstrate that isoprene is the dominant HCHO precursor even in these cases (Sect. 5). 

160 Thus, chemistry (and not co-variance of NO x and anthropogenic VOC) drives the observed NO x 

161 dependence of HCHO abundance. 

162 

163 3. Linking Observed and Emitted Isoprene 

164 The isoprene photochemical cascade is a multi-step process. Isoprene oxidation is initiated by reaction 

165 with the hydroxyl radical (OH), ozone or the nitrate radical (NO3). In the Southeast U.S., typical daytime 

166 levels for OH, ozone and NO 3 are 4 x 10 6 cm' 3 , 50 ppbv and 0.1 pptv, respectively (OH and N0 3 are 

167 estimated from median box model output, see Sect. 5). The corresponding isoprene lifetimes at 298K 

168 are 0.7 h, 17 h and 160 h, respectively. Thus, reaction with OH typically constitutes 95% or more of the 

169 total daytime isoprene sink in this environment. Addition of OH and reaction with 0 2 generates one of 

170 several isoprene hydroxyperoxy radicals (ISOP0 2 ). ISOP0 2 isomers interconvert rapidly due to reversible 

171 0 2 addition (Peeters et al., 2009) but are eventually destroyed via reaction with NO, hydroperoxy radical 

172 (H0 2 ), other organic peroxy radicals (R0 2 ) or isomerization. Most branches have the potential to 

173 produce HCHO, with varying yields. The laboratory-derived first-generation HCHO yield from the NO 

174 pathway is ~0.6 (Atkinson and Arey, 2003), though this value may be less representative of the real 

175 atmosphere due to the very high isoprene concentrations (and very short R0 2 lifetimes) in early 

176 chamber experiments. The first-generation yield from the H0 2 pathway is ~0.06 (Liu et al., 2013). 
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177 Isomerization chemistry is less well understood; the 1,5-H-shift is believed to produce HCHO with a unity 

178 yield, while the much faster 1,6-H-shift should not produce any HCHO (da Silva et al., 2010; Fuchs et al., 

179 2013; Peeters et al., 2014; Peeters and Muller, 2010; Peeters et al., 2009). Regardless of the specific 

180 pathway, MVK or MACR are always co-produced with HCHO in the first generation. HCHO is also 

181 generated in subsequent chemistry, but on a longer timescale and from a much larger suite of 

182 precursors. For example, the OH lifetimes of MACR and MVK are respectively 3.5 and 5 times longer 

183 than that of isoprene. 

184 Boundary layer composition reflects a mixture of emissions with various degrees of 

185 photochemical processing. To isolate the impact of "fresh" isoprene emissions, we exploit the relatively 

186 simple chemistry of MVK and MACR, which are produced via isoprene oxidation and lost primarily via 

187 reaction with OH. 


188 

ISOP + OH YmacrMACR + YmvkMVK 

ki = 2.7xl0 11 e 390/T 

(Rl) 

189 

MACR + OH products 

k 2 = 8.0xl0 12 e 380/T 

(R2) 

190 

MVK + OH products 

k 3 = 2.6xl0 12 e 610/T 

(R3) 


191 Rate constants (k) are taken from the IUPAC database (Atkinson et al., 2006). These reactions form the 

192 basis for a photochemical clock of isoprene oxidation (de Gouw et al., 2005; Roberts et al., 2006; Stroud 

193 et al., 2001). Integration of the kinetic equations for this system shows that the product/parent ratios 

194 are a function of the rate constants, yield (y), reaction time (f) and the mean OH concentration averaged 

195 over reaction time. In the case of MACR, for example: 

196 [MACR] = VMACRk 1 ^ c (ki-k^\OH\t\ 

197 An analogous expression holds for MVK. As noted by Stroud et al. (2001), this "sequential reaction 

198 model" is purely chemical and does not account for the effects of mixing and transport. Indeed, this 

199 analysis relates daughter/parent ratios to an "average" photochemical age, when in fact there is a broad 

200 distribution of ages in any mixed air mass. We also implicitly assume that direct emissions (Fares et al., 

201 2015) and deposition (Karl et al., 2010) of MVK and MACR do not significantly influence the budget of 

202 these compounds. 

203 Two potential issues arise when applying this model to the real atmosphere. First, the yields of 

204 MVK and MACR are dependent on ISOP0 2 branching and are thus a non-linear function of NO x . Previous 
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205 applications of this method (de Gouw et al., 2005; Roberts et al., 2006; Stroud et al., 2001) have 

206 assumed lab-derived high-NO x yields of 0.33 and 0.23 for MVK and MACR, respectively (Atkinson and 

207 Arey, 2003), but this may not be appropriate in the present case; furthermore, these yields are not fully 

208 consistent with current chemical mechanisms (Fig. S4). We explicitly examine the effects of NO x -varying 

209 yields below using yield curves derived from box model simulations (see SI for details). Second, the 

210 photochemical age (f) implied by any observed daughter/parent ratio depends on the concentration of 

211 OH, which was not measured and varies as an air mass ages. Rather than assume a single "typical" value 

212 for OH, we express photochemical age in terms of "exposure," defined here as the product of OH 

213 concentration and reaction time integrated over the photochemical lifetime of an air mass. 

214 Figure 2 compares the observed relationship of MVK/isoprene and MACR/isoprene ratios 

215 against theoretical trends predicted by the sequential reaction model. Theoretical ratios are calculated 

216 at fixed exposures of 2, 4, 8, 12 and 16x 10 6 OH cm' 3 h using two sets of yields: high NO (NO = 1000 pptv, 

217 y M vK = 0.41, yMACR = 0.28) and low NO (NO = 50 pptv, y mvk = 0.21, y M ACR = 0.19). Observed ratios of 

218 MVK/isoprene versus MACR/isoprene exhibit a tight linear correlation (Fig. 2). Higher ratios are often 

219 associated with higher NO x , likely reflecting enhanced OH and higher HCHO yields in these air masses. 

220 Far downwind from isoprene and NO x source regions, we would expect to see higher MVK/isoprene and 

221 MACR/isoprene ratios associated with lower NO x due to removal of the latter. The theoretical slope 

222 agrees well with observations, indicating exposures of 1 - 16 x 10 6 OH cm' 3 h. For a typical daytime OH 

223 concentration of 4 x 10 6 cm' 3 , this corresponds to processing times of 15 minutes to 4 hours. 

224 The ratio of y M vK to Ymacr dictates the location of the theoretical line and thus the 

225 correspondence between daughter/parent ratios and exposure. For example, a MACR/isoprene ratio of 

226 1 would be consistent with an exposure of 4.9 x 10 6 OH cm' 3 h at high-NO x conditions (NO = 1000 pptv) 

227 versus 6.1 x 10 6 OH cm 3 h at low-NO x (NO = 50 pptv). Thus, for any given daughter/parent ratio, a higher 

228 assumed yield gives a smaller derived exposure. Observations in Fig. 2 fall above the high-NO x 

229 theoretical relationship. As discussed in the SI, however, iWAS MVK measurement may contain a 

230 positive artifact on the order of 50%. This potential systematic error (thick black line in Fig. 2) overlaps 

231 both the high and low-NO x theoretical relationships. Given the wide range of conditions sampled, it is 

232 most appropriate to use a NO x -dependent yield for MVK and MACR. For this purpose, model-derived 

233 yields (Fig. S4 and SI) are interpolated to observed NO mixing ratios. 
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234 We can effectively reverse this photochemical clock to derive a proxy for the total isoprene 

235 emissions that had been released into the sample air masses (de Gouw et al., 2005). First, we calculate 

236 OH exposures from observed daughter/parent ratios by inverting Eqn. (1). To perform this calculation 

237 with PTR-MS data (which has far greater coverage), we partition the measured sum between MVK and 

238 MACR using MVK/MACR ratios from box model calculations (Sect. 5). Modeled MVK/MACR ratios (with 

239 an output interval of 1 minute) are linearly interpolated to the 14-second observational time base. The 

240 MVK/MACR ratio does not vary dramatically (mean ± Ict: 1.3 ± 0.2), and using a constant ratio instead 

241 alters results by less than 4%. Calculated exposures range from 0.5 to 18 x 10 s OH crrf 3 h (Fig. S5A). 

242 Exposures derived from MACR are 6% higher than those from MVK on average, and we use the mean of 

243 these two values. Next, an "initial" isoprene mixing ratio, ISOPo, is estimated via reverse integration of 

244 isoprene's first-order loss rate: 

245 [ ISOP] 0 = [IS0P]e k ^ OH]t (2) 

246 ISOPo represents the amount of isoprene that an air parcel would have to start with to generate the 

247 amount of isoprene, MVK and MACR observed. Thus, it is an observationally-constrained surrogate for 

248 isoprene emission strength (modulated to some degree by boundary layer height, as it is a volume- 

249 based quantity). ISOPo mixing ratios are typically 2-10 times higher than observed isoprene (Fig. S5B). 

250 

251 4. The Yield of HCHO from Isoprene 

252 The definition of "yield" can vary with context and requires careful consideration when quantifying the 

253 isoprene-HCHO relationship. In a mechanistic sense, the "first generation yield" refers to the amount of 

254 HCHO produced per unit isoprene consumed in the first stage of oxidation. This is analogous to the 

255 yields of MVK and MACR used in the above calculation of initial isoprene. The model-derived first- 

256 generation HCHO yield from isoprene varies by more than a factor of 2 over the range of chemical 

257 environments encountered during SENEX (Fig. S4). An alternative definition is that of the "total yield" 

258 (sometimes referred to as the "molar yield," e.g. Millet et al. (2006)), a time-dependent quantity that 

259 describes the total amount of HCHO produced over multiple generations of oxidation. The total yield is 

260 typically derived from model simulations and used to relate satellite HCHO column observations to 

261 isoprene emissions (Marais et al., 2012; Millet et al., 2006). Early studies acknowledged the NO x 

262 dependence of the total yield (Millet et al., 2006; Palmer et al., 2003) and more recent work has 
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263 attempted to account for the dependence using NO 2 column observations (Marais et al., 2012). Here, 

264 we define the "prompt yield" as the change in observed HCHO per unit change in ISOP 0 (y P = 

265 AHCHO/AISOPo). This is not the same as the first-generation yield, since y p can include HCHO production 

266 and loss over several hours (depending on the photochemical exposure of an air mass). Nor is it the 

267 same as the total yield, which inherently does not account for HCHO loss as an air mass ages. The 

268 prompt yield is effectively a quantity that relates isoprene emission strength to observed HCHO 

269 abundance. As we will demonstrate, y p is well-suited for segregating the various drivers of HCHO and for 

270 benchmarking model performance. 

271 Figure 3A shows the relationship between calculated ISOP 0 and observed HCHO. The overall 

272 correlation is linear with a striking NO x gradient. To quantify this NO x dependence, we sort the data by 

273 log(NO x ), group it into 20 bins such that each bin contains the same number of points (N = 416), and 

274 perform a major-axis linear fit of HCHO versus ISOP 0 for each bin. Individual fits give r 2 values of 0.6-0. 8, 

275 except for the highest NO x bin (r 2 = 0.48) that contains some heavily-polluted air masses, such as 

276 downwind from power plants. Results are independent of the number of bins chosen or time resolution 

277 (e.g., 1-second versus 1-minute data). 

278 The HCHO-ISOP 0 slope (Fig. 3B) represents the prompt yield. This yield varies by a factor of 3 

279 over the range of observed NO x , from 0.3 for NO x mixing ratios of a few hundred pptv to 0.9 at NO x > 1 

280 ppbv. At low NO x , y p is comparable to the MCM-predicted direct first-generation yield of HCHO (0. 3-0.4 

281 at NO = 10-40 pptv, Fig. S4), while at high NO x it is somewhat higher than the predicted first-generation 

282 yield (0.74 at NO = 1000 pptv). This likely reflects the inclusion of more than one generation of HCHO 

283 production at higher NO x , where oxidation is more rapid (median exposures increase by 38% over the 

284 range of observed NO x values). Most of this portion of the HCHO budget, however, stems from first- 

285 generation production. 

286 The intercept (Fig. 3C) represents the abundance of "background" HCHO. This portion of the 

287 HCHO budget stems mainly from air that either has not encountered strong isoprene emissions or is so 

288 aged that most of the isoprene has reacted away and can no longer be linked to a specific source region. 

289 Some of this background may also stem from oxidation of long-lived primary emissions like methane or 

290 methanol. Box model calculations (Sect. 5) indicate average HCHO budget contributions of 0.3 ± 0.2 

291 ppbv and 0.2 ± 0.1 ppbv from methane and methanol, respectively. Background HCHO also exhibits a 

292 marked NO x dependence, increasing from 1.6 to 3.3 ppbv over the observed NO x range. As with y p , we 
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293 expect such behavior since NO x regulates the fate of all organic peroxy radicals (see Sect. 6). Assuming a 

294 1 km mixed layer depth (Wagner et al., 2015), the corresponding HCHO column density for this 

295 background is 4 - 8 x 10 15 cm' 2 . This is comparable to the background reported by previous 

296 investigations of satellite-derived HCHO columns (Barkley et al., 2013; Millet et al., 2006). None of these 

297 studies explicitly account for the NO x dependence of the background, though it can represent a 

298 substantial fraction of the total HCHO column - maximum summertime HCHO columns over the 

299 southeast U.S. are ~25 x 10 15 cm' 2 (Millet et al., 2008). Given the strong NO x dependence of both y p and 

300 background HCHO, grouping HCHO column observations by NO x (e.g. using simultaneous observations of 

301 N0 2 columns (Marais et al., 2012) or model-derived NO x ) and performing an analysis similar to that 

302 described here should provide a robust means of accounting for these influences. 

303 

304 5. Model Evaluation 

305 To illustrate the utility of this analysis, we compare the observed HCHO-ISOP 0 relationship to results 

306 from a global chemical-transport model and a 0-D box model. Goals? 

307 The GFDL AM3 model is an atmospheric general circulation model with interactive chemistry 

308 (Donner et al., 2011), including recent updates to the representation of isoprene degradation (Mao et 

309 al., 2013; Naik et al., 2013). Model simulations were carried out at 50 x 50 km 2 resolution with 

310 horizontal winds nudged to NCEP GFS analyses and sampled along the SENEX flight tracks at a time 

311 resolution of 1 minute. Further details are available elsewhere (Li et al., 2015). 

312 The University of Washington Chemical Box Model (UWCM v2.2) is a versatile 0-dimensional 

313 framework for simulating various chemical systems, including lab chamber experiments (Wolfe et al., 

314 2012) and observations from ground (Kim et al., 2015; Kim et al., 2013; Wolfe et al., 2014) and airborne 

315 (Marvin et al., 2015) platforms. Multiple chemical mechanisms are available within UWCM; here we use 

316 the latest version of the Master Chemical Mechanism (MCM v3.3, Jenkin et al. (2015)). UWCM is 

317 constrained with 1-minute average observations of isoprene, N0 2 , ozone, CO, PAN, methane, methanol 

318 and meteorology and assumes clear-sky conditions for photolysis frequencies. The chemical system is 

319 integrated forward in time to diel steady state (total integration time of 3 days) for each set of 

320 measurements. This setup inherently assumes that the atmosphere is in chemical steady state - that is, 

321 that production and loss of HCHO, MVK, MACR and other species are roughly balanced. This assumption 

322 is rarely strictly true and may fail for highly-aged air masses (where isoprene is depleted) or close to 
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323 strong local emissions. Nonetheless, it is a fair approximation for the daytime well-mixed boundary layer 

324 observations that prevailed during SENEX. Monoterpenes and anthropogenic VOC are excluded from the 

325 simulation since observations of these species (from the iWAS) are relatively sparse. Separate sensitivity 

326 simulations utilizing the iWAS data suggest that observed monoterpenes and anthropogenic VOC (a 

327 subset of alkanes, alkenes and aromatics) increase modeled HCHO by 1 ± 2 % and 2 ± 3 %, respectively. 

328 A more detailed evaluation of box model performance is forthcoming (Marvin et al., 2015). 

329 Output from both models is filtered for daytime, boundary-layer, non-biomass burning points 

330 using the same criteria as that for observations (Sect. 2). Both models adequately reproduce observed 

331 HCHO mixing ratios (Fig. S6). We perform the same analyses as described above to derive model y p and 

332 background HCHO. Because of the reduced time resolution, we group results into 10 NO x bins, instead of 

333 20, before fitting. For AM3, this results in 172 points per bin and typical r 2 values of 0.5 - 0.8. For 

334 UWCM, there are 157 points per bin and all r 2 values are > 0.9. 

335 Both AM3 and UWCM reproduce the observed NO x dependence of the prompt yield (Fig. 4A). 

336 AM3 agrees well with observations in both magnitude and trend, though with some scatter at mid-NO x 

337 levels. UWCM tends be slightly high throughout the whole NO x range, which may reflect an over- 

338 estimation of first-generation HCHO production due to holding isoprene constant throughout the model 

339 step and/or assuming diel steady state. Regardless, these results suggest that both models provide 

340 excellent representation of early generation isoprene oxidation across NO x regimes. 

341 Background HCHO mixing ratios are under-predicted by 0.5 - 1 ppbv by both models (Fig. 4B). 

342 The range of under-prediction is consistent with the offsets between observed and modeled total HCHO 

343 abundances (Fig. S6 fit x-intercepts: 0.3 ppbv (AM3) and 0.9 ppbv (UWCM)). It is possible that both 

344 models are missing some HCHO precursors (e.g. from multi-generation isoprene oxidation or other VOC 

345 not related to isoprene). This is especially plausible for the UWCM simulation, which only includes 

346 isoprene, methane and methanol as primary VOC and does not account for horizontal transport. Under- 

347 estimated OH concentrations might also explain part of this discrepancy, though we cannot easily 

348 evaluate this possibility. AM3 performs somewhat better than UWCM in terms of overall magnitude but 

349 exhibits a less clear NO x trend, which may reflect dilution over fairly large grid scales (note that the 

350 range of binned NO x values is smaller for AM3 than both observations and the UWCM). This result again 

351 highlights the need to consider this background before using a model to interpret observed HCHO 

352 columns that effectively integrate HCHO sources over space and time. 
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353 


354 6. Mechanistic Drivers of the NO x - HCHO Relationship 

355 Despite the complexity of gas-phase organic chemistry, the impact of NO x on HCHO production 

356 essentially reduces to two factors: radical cycling and R0 2 branching. Increasing NO enhances the 

357 conversion of H0 2 to OH (R4) and thus accelerates VOC oxidation (R5) and HCHO loss. Subsequent 

358 production of HCHO depends on the structure and fate of R0 2 intermediates, which can react with NO, 

359 H0 2 , other R0 2 , or isomerize (R6). 

360 N0 + H0 2 ^ N0 2 + 0H (R4) 

361 V0C + 0H^R0 2 (R5) 

362 R0 2 + (NO, H0 2 , R0 2 , isomerization) -> aHCHO (R6) 

363 Here, a represents a bulk branching ratio for HCHO production weighted over all R0 2 reactions. The R0 2 

364 lifetime is typically less than 100 s during the day, so (R5) is the rate-limiting step in HCHO formation. 

365 The HCHO production rate is then equal to the product of the total R0 2 production rate and the bulk 

366 branching ratio. 

367 To disentangle these factors, we extract chemical rates from the diel steady-state UWCM 

368 simulations discussed in Sect. 5. Figure 5 shows the gross production rates of total peroxy radicals and 

369 HCHO as a function of NO x . Consistent with our earlier discussion, total HCHO production increases by 

370 more than a factor of 3 from low to high NO x . In contrast, R0 2 production is effectively constant within 

371 model variability. Closer investigation (results not shown) reveals that a factor of 3 - 4 increase in OH 

372 concentrations between low and high NO x is more than offset by a similar reduction in isoprene. The 

373 ratio of HCHO to R0 2 production rates gives an estimate for a, which increases from 0.14 to 0.39 across 

374 this NO x range (Fig. 5). Though the total R0 2 production rate includes reactions that do not make HCHO, 

375 a is still a useful metric for the relationship between HCHO production and overall VOC oxidation. Based 

376 on this analysis, we conclude that changes in R0 2 branching are the dominant factor driving the NO x 

377 dependence of HCHO production and abundance. 

378 Increased OH also reduces the lifetime of HCHO, which may affect the HCHO budget if this 

379 reaction becomes competitive with photolysis. UWCM predicts an average HCHO photolysis lifetime of 4 

380 hours and OH reaction lifetimes that range from 3 hours at high NO x to 12 hours at low NO x . Thus, 
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381 photolysis is typically the dominant loss process and the scaling of HCHO lifetime with OH is typically 

382 weak. As a result, the net chemical tendency of HCHO (production minus loss, not shown) is positive and 

383 increasing throughout the range of model NO x conditions. Faster loss due to reaction with OH therefore 

384 only slightly dampens the enhancement in HCHO production. 

385 

386 7. Conclusions 

387 Using SENEX aircraft observations, we have quantified the NO x dependence of the relationship between 

388 isoprene emission strength and HCHO mixing ratios. Simultaneous measurements of isoprene, MVK and 

389 MACR define a photochemical clock for isoprene oxidation, allowing separation of prompt HCHO 

390 production (which retains the isoprene source signature) and background HCHO from late-generation 

391 isoprene oxidation products, methane and other long-lived VOC. The prompt HCHO yield increases by a 

392 factor of 3 (0.3 to 0.9 mol/mol) and the average background HCHO mixing ratio more than doubles (1.6 

393 to 3.3 ppbv) over the range of NO x values encountered in the southeast U.S. This analytical method is 

394 applied to evaluate the performance of a global chemical transport model and a 0-D box model. Both 

395 models accurately reproduce the observed NO x trend of the prompt HCHO yield, indicating that both 

396 chemical mechanisms accurately capture early-stage isoprene oxidation. On the other hand, both 

397 models also under-predict background HCHO abundance by 0.5 - 1 ppbv, which may be a significant 

398 fraction of total HCHO in some cases. This may suggest insufficient build-up of isoprene-derived long- 

399 lived precursors in the models, missing VOC not related to isoprene, or insufficient OH. Box model 

400 results also provide insight into the mechanistic drivers of the observed NO x trends. We find that 

401 increasing NO x does not significantly affect total R0 2 production due to the cancelling effects of higher 

402 OH and lower VOC, and thus the positive correlation between NO x and HCHO primarily reflects the 

403 changing fate of R0 2 radicals. 

404 To our knowledge, there are no direct laboratory measurements of HCHO yields from low-NO x 

405 isoprene chemistry; thus, the results presented here constitute the first measurement-constrained 

406 evaluation of the isoprene-HCHO link across NO x regimes. The AM3 and MCMv3.3 mechanisms differ 

407 substantially (the former is highly condensed while the latter is explicit), but both contain recent 

408 updates to isoprene degradation. We expect that other mechanisms will also perform well if they 

409 accurately reflect our current best understanding. The observations presented here do not include the 

410 extremely-low NO x regime (NO x < 0.1 ppbv) typical of remote regions like the Amazon and equatorial 
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411 Africa. In such pristine regions, smearing of HCHO production is expected to be more severe (Barkley et 

412 al., 2013), and total HCHO production may be significantly lower if the R0 2 fate favors functionalization 

413 over fragmentation (e.g. isomerization). More work is needed to map out this area of the urban-rural 

414 spectrum. It may also be possible to apply the methods developed here to evaluate the chemistry of 

415 glyoxal, another key tracer of VOC oxidation that is also amenable to orbital observations (Kaiser et al., 

416 2015b; Li et al., 2015) and is believed to be an important precursor for SOA (McNeill et al., 2012). 

417 These results also carry implications for top-down isoprene emission estimates. Uncertainties in 

418 low-NOx chemistry are often cited as the largest source of potential error in derived emissions (Marais et 

419 al., 2012; Palmer et al., 2006). Based on our analysis, current mechanisms appear to capture low-NO x 

420 production of HCHO, MVK and MACR, thus such errors are likely less severe than commonly asserted. 

421 Recent work has acknowledged the impact of NO x on the prompt yield of HCHO from isoprene (Marais 

422 et al., 2012). We advocate considering the NO x dependence of background HCHO as well, since this can 

423 constitute a significant fraction of the total HCHO column. For scale, the derived background HCHO 

424 mixing ratio of 1.6 - 3.3 ppbv is 37 - 77% of the campaign-mean observed HCHO mixing ratio of 4.3 

425 ppbv. Forthcoming geostationary observations will deliver sufficient resolution to delineate local 

426 gradients in chemical regime, and smearing and background HCHO production will become problematic 

427 even in high-NO x regions. Indeed, even current-generation orbital instruments are capable of resolving 

428 urban-rural gradients in HCHO columns (Boeke et al., 2011). When applying advanced statistical 

429 techniques like inversion, model results will only be as accurate as the chemical mechanisms driving 

430 them. Continued field observations are crucial for providing confidence in our ability to link HCHO to its 

431 sources. 

432 
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677 Figure 1 . Co-variation of isoprene, NO and HCHO mixing ratios in the summertime southeast U.S. Data 

678 are limited to daytime boundary layer observations. Histograms show the sampled NO and isoprene 

679 distributions. 
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Figure 2. A photochemical clock of isoprene oxidation defined by the progression of daughter/parent 
ratios. Solid circles show the observed ratios calculated from iWAS observations, colored by NO x . 
Blue/purple symbols, dashed lines, and text indicate the theoretical exposures (the product of OH 
concentration and time) corresponding to any given daughter/parent relationship. Theoretical values 
are calculated at 298K using MVK and MACR yields for NO values of 50 pptv (triangles) and 1000 pptv 
(squares). The thick black line denotes the systematic error due to a potential 50% positive artifact in 
MVK observations (see SI). 
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Figure 3. (A) NO x modulates the relationship between observed HCHO and calculated initial isoprene 
mixing ratios. Symbols denote all 1-second data points. Dashed lines illustrate representative major-axis 
fits of NOx-grouped subsets at mean NO x values of 170, 380 and 810 pptv (see text for details of fitting 
procedure). The slope (B) and intercept (C) of these fits are the prompt HCHO yield and background 
HCHO mixing ratio, respectively. Error bars in (B) and (C) are 3a fitting uncertainties. 
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Figure 4. Comparison of observed and model-derived relationship between HCHO and initial isoprene 
versus NO x . Slopes (A) and intercepts (B) are calculated as described in the text. The observed values 
(blue line with shading) are the same as those shown in Figs. 3B-C. Symbols represent fit results for the 
global AM3 model (red circles) and the 0-D UWCM box model (black diamonds). Error bars denote 3a 
fitting uncertainties. 
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704 Figure 5. NO x dependence of chemical properties related to HCHO production, extracted from the 

705 UWCM simulation of SENEX observations. Production rates for HCHO (blue) and total R0 2 (orange) are 

706 averaged over NO x using 10 bins with equal numbers of points. Solid lines show the mean, shading is la 

707 variability. Note that R0 2 production is scaled down by a factor of 10. The ratio of HCHO to R0 2 

708 production gives the bulk HCHO branching ratio (dashed line). 
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